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We report on a first search for resonant pair production of neutral long-lived particles (NLLP) 
which each decay to a bb pair, using 3.6 fb _1 of data recorded with the DO detector at the Fermilab 
Tevatron collider. We search for pairs of displaced vertices in the tracking detector at radii in the 
range 1.6-20 cm from the beam axis. No significant excess is observed above background, and upper 
limits are set on the production rate in a hidden- valley benchmark model for a range of Higgs boson 
masses and NLLP masses and lifetimes. 

PACS numbers: 12.60.Fr, 14.80.Cp 



A class of hidden-valley (HV) models [l[ predicts a 
new, confining gauge group that is weakly coupled to the 
standard model (SM), leading to the production of HV 
particles (v-particles). The details of v-particle decay de- 



pend on the specific model, but the HV quarks always 
hadronize due to confinement producing "v-hadrons" 
that can be long-lived. One particular model used as 
a benchmark for this search is the SM Higgs boson (H) 
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mixing with a HV Higgs boson that gives mass to v- 
particles. The SM Higgs boson could then decay di- 
rectly to v-hadrons through this mixing with a substan- 
tial branching fraction Q. These v-hadrons may couple 
preferentially to heavy SM particles, such as b quarks, 
due to helicity suppression. The result is a striking ex- 
perimental signature of highly displaced secondary ver- 
tices (SV) with a large number of attached tracks from 
the b quark decays. Direct searches at the CERN LEP 
collider have excluded a Higgs boson decaying to bb or 
rf with M H < 114.4 GeV at the 95% C.L. [|. But if the 
Higgs boson dominantly decays to long-lived v-particles 
which then decay inside the detector to bb, only the most 
general LEP limit is relevant, Mjj > 81 GeV, for any 
Higgs boson radiating off a Z boson Cosmological 
constraints require that one of the light v-hadrons have 
a lifetime <C 1 second to be consistent with models of big 
bang nucleosynthesis 

In this Letter, we present the first search for pair- 
produced neutral long-lived particles (NLLP), each de- 
caying to a b quark pair, using the DO detector [j| at the 
Fermilab Tevatron pp collider. The b quarks are required 
in order to provide a high transverse momentum (pt) 
muon for triggering with high efficiency The data were 
collected from April 2002 to August 2008 and correspond 
to an integrated luminosity of 3.6 fb _1 at ^/s = 1.96 TeV. 
The DO central tracking detector comprises a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT), 
both located within a 2 T superconducting solenoidal 
magnet. The SMT, extending from a radius of «2 cm 
to «10 cm, has a six-barrel longitudinal structure, each 
with a set of four layers arranged axially around the 
beam pipe, and interspersed with 16 radial disks. The 
CFT, extending from a radius of «20 cm to ~50 cm, has 
eight thin coaxial barrels, each supporting two doublets 
of overlapping scintillating fibers. Secondary vertices are 
reconstructed by combining charged particle tracks found 
in the tracking detector, which effectively limits the anal- 
ysis to NLLP decays occurring within a maximum radius 
of 20 cm, well within the tracker volume. We also exclude 
vertex radii less than 1.6 cm since the background from 
heavy flavor production is large in that region. Known 
sources of SVs other than heavy-flavor include decays in- 
flight of light particles, inelastic interactions of particles 
with nuclei of detector material, and photon conversions. 
Vertices may also be mimicked by pattern recognition 
errors. 

PYTHIA [6j is used to simulate signal and background 
events, which are then passed through a full GEANT3- 
based 0] DO detector simulation and the same recon- 
struction as for collider data. For signal, the SM gg^H 
process is generated, the Higgs boson is forced to decay 
to a pair of long-lived A bosons (a heavy, neutral scalar, 
representing a v-hadron), and each A boson is forced to 
decay to a pair of b quarks. The Higgs boson mass (M#) 
is varied from 90 to 200 GeV, the v-hadron mass (rriHv) 
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FIG. 1: Material map in the plane transverse to the beam- 
line, generated using SVs with three attached tracks, in events 
passing the initial selection. The structure of the silicon de- 
tector and supports are clearly seen. 

from 15 to 40 GeV and the average v-hadron proper de- 
cay length (Ld = or) from 2.5 cm to 10 cm. For back- 
ground, inclusive pp multijet events are generated. Ap- 
proximately one hundred thousand Monte Carlo (MC) 
events for each signal sample and ten million events of 
multijet background are generated and are overlaid with 
data to simulate detector noise and pile-up effects from 
additional pp interactions. 

At least two jets with a cone radius of 0.5 § are re- 
quired, each with pt > 10 GeV. And at least one muon 
is required with pr > 4 GeV, matched within A1Z < 0.7 
to one of the jets, where A1Z — y/ (Ac/)) 2 + (Ar/) 2 with 
(f> being the azimuthal angle and r\ the pseudorapidity. 
The muon requirement is more efficient for signal than 
background due to the presence of a b — > fi or b — ► c — > fi 
decay from at least one of the four b quarks, and is also 
required for an accurate measurement of the trigger ef- 
ficiency. Primary vertices (PVs) are reconstructed by 
clustering tracks and correspond to pp interaction lo- 
cations. To ensure good SV reconstruction, we further 
require fewer than four PVs be reconstructed and that 
the selected PV with the largest JV logp^, summed over 
all vertex tracks i, be located within \z\ < 35 cm and 
r < 1 cm, where x and y are the horizontal and vertical 
components of the distance r with respect to the beam 
axis, and z is the distance along the beam axis from the 
center of the detector. An initial selection requires that 
each event has at least one SV with 2D decay length from 
the PV in the plane transverse to the beam (L x d v ) larger 
than 1 cm and decay length significance (decay length 
divided by its uncertainty) greater than five. The mo- 
mentum of the SV, reconstructed from the vectorial sum 
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of the momenta of its associated tracks, must point away 
from the PV to reduce combinatoric background. SVs are 
reconstructed using a track selection so as to efficiently 
combine the b and b decay products of each v-hadron into 
a single SV. Approximately 50 million data events satisfy 
these requirements, dominated by dijet and heavy-flavor 
production. 

To maximize the discovery potential of this analysis we 
use an OR of all triggers. The most frequently fired trig- 
gers that make up the datasct passing the initial selection 
involve a muon and jet at the first trigger level and re- 
finements of these objects at higher levels. The overall 
trigger OR efficiency is estimated by first measuring the 
efficiency for a single trigger per data collection period 
using known muon and jet trigger efficiencies. Then the 
number of events fired by that single trigger is compared 
to the total number of data events passing the OR of 
all triggers, as a function of sensitive variables, such as 
muon pt, jet pr, jet angles, etc. No significant depen- 
dence is found, except on jet px, thus the overall trigger 
OR efficiency is modeled as a function of jet pr- 

Further selections are optimized by maximizing the ex- 
pected signal significance (S/y/S + B), where B and S 
are the number of MC background and signal events, 
respectively. The heavy-flavor background, mainly b 
hadrons with cr«0.3 cm, produces a very large number 
of SVs, but their number decreases exponentially as the 
radial distance of the SV from the PV increases. SVs are 
required to have L^ v >1.6 cm. We expect signal events 
to preferentially produce SVs with a large number of at- 
tached tracks, therefore we require SVs with track mul- 
tiplicity of at least four. Interactions of primary collision 
particles (tt, protons, etc.) with detector material, such 
as silicon sensors, cables, etc., are the major source of 
background. In order to quantify the material regions, 
we construct a map of SV density in data, using SV with 
track multiplicity of three, in the xy (see Fig. [1]) and 
rz projections. SVs that occur in regions of high SV 
density are then removed. After this "pre-selection" is 
performed, the multijet background MC sample is nor- 
malized to the data (see Table QJ. Finally, at least two 
SVs are required in each event, and they are required to 
have ATZ(SV1, SV2)>0.5, to prevent cases where a single 
true vertex is mis-reconstructed as two nearby separate 
vertices. No events in data have more than two SVs. 

Two more variables are used to select the signal: SV 
invariant mass and SV collinearity. The invariant mass 
is reconstructed from the four-momenta of the outgoing 
tracks attached to a SV, assuming the pion mass for all 
particles. Collinearity is defined as the cosine of the angle 
between the vector sum of the momenta of the attached 
tracks and the direction to the SV from the PV. De- 
pending on the signal point, one of these two variables is 
used to perform the final separation of signal and back- 
ground. The quality of the background model is of pri- 
mary importance so we develop a method of tuning the 



TABLE I: Summary of event selections, showing the remain- 
ing background, data, and signal events (for Mjf=120 GeV 
and Ld=5 cm) after each selection. Background is normalized 
to data after pre-selection. 
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FIG. 2: The minimum mass of the two SVs, for data, 
background MC, and signal MC with M H =120 GeV, 
iriHV=15 GeV, and Ld=5 cm. The hatched region shows 
the uncertainty on the background MC. 
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FIG. 3: The maximum collinearity of the two SVs, for 
data, background MC, and signal MC with A/h = 120 GeV, 
m_frv r =40 GeV, and Ld~5 cm. The hatched region shows the 
uncertainty on the background MC. 
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multijet background simulation of the SV invariant mass 
and SV collinearity distributions to data. The events af- 
ter pre-selection are divided into two distinct sets: the 
first contains events with only one SV (1SV), whereas 
the second contains events with at least two SVs (2SV). 
Since the signal content of the 1SV set is expected to be 
<0.1% for any of the signal MC points studied, we use the 
1SV set to compare the data to the multijet background 
MC and perform corrections to the MC. Gaussian smear- 
ing functions are applied to fractions of the background 
MC events for the SV invariant mass and SV collinearity 
to model the tails of the distributions better. The SV 
invariant mass is smeared using a width of 12 GeV in 
about 1% of events, and the SV collinearity is smeared 
with a width of 0.15 in about 1.5% of events. The same 
smearing is then also applied to all MC signal samples. 
For rriHv < 20 GeV, a requirement on the minimum SV 
mass in an event >2.5 GeV is most effective (Fig. [2]). 
For heavier v-hadrons, we take advantage of the SVs de- 
cay products being more widely spread in angle, which 
is better measured than the invariant mass. Requiring 
the maximum SV collinearity in an event to be <0.9937 
maximizes the expected significance (Fig. [3]). 

The uncertainty on the signal acceptance is dominated 
by the modeling of trigger efficiency and is (13-17)%. 
The uncertainty on the background due to the difference 
in track reconstruction efficiency between MC and data 
is estimated by using two methods of normalization and 
found to be 28%. We estimate the effect of smearing the 
MC samples by performing the entire analysis without 
smearing. For the requirements applied to the SV mass 
or collinearity, smearing results in a difference of up to 
18% on the multijet background yield and a negligible 
difference on the signal acceptance. Smearing also has 
no effect on the optimized requirement values. To esti- 
mate the uncertainty from requiring a small SV density, 
we compare the difference in the number of remaining 
events between multijet background and data before and 
after making the density requirement, and find agree- 
ment within (8-15)%. The uncertainty on the integrated 
luminosity is 6.1% 

The final results after all selections are sum- 
marized in Table HU No significant excess 
is observed, so 95% C.L. limits are set on 
a(H+X)xBR{H->HVHV)xBR 2 {HV->bb) using a 
modified frequentist method [HI ], which includes all sys- 
tematic uncertainties on signal acceptance, background, 
and luminosity. Depending on the signal parameters, 
Higgs boson production about 1-10 times the SM cross 
section is excluded, if the Higgs boson always decays to 
a pair of long-lived v-hadrons decaying only to bb (see 
Fig. [4j. These results also provide the first constraints 
on pair-produced NLLPs decaying to b jets in the radial 
range of 1.6-20 cm at a hadron collider. 

We thank the staffs at Fermilab and collaborating 
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100% for~BR{H-*HVHV) and BR{HV~*bb). (color online) 
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TABLE II: Results for each simulated signal: the numbers of background, signal, and data events after all selections, overall 
signal efficiency, SM Higgs production rate, and observed and expected 95% C.L. upper limits on the signal cross section. 
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